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Abstract 

Sufficiently small Schwarzschild black holes in global AdSsxS^ are Gregory-Laflamme 
unstable. We construct new families of black hole solutions that bifurcate from the onset of 
this instability and break the full SO(6) symmetry group of the S^ down to SO(5). These 
new “lumpy” solutions are labelled by the harmonics We find evidence that the ^ = 1 
branch never dominates the microcanonical/canonical ensembles and connects through a 
topology-changing merger to a localised black hole solution with S^ topology. We argue that 
these S^ black holes should become the dominant phase in the microcanonical ensemble for 
small enough energies, and that the transition to Schwarzschild black holes is first order. 
Furthermore, we find two branches of solutions with £ = 2. We expect one of these branches 
to connect to a solution containing two localised black holes, while the other branch connects 
to a black hole solution with horizon topology S^ X S^ which we call a “black belt”. 
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1 Introduction 

Gauge-gravity duality provides a unique arena to study quantum gravity in its most extreme 
regimes. Its best understood formulation stipulates an equivalence between ten dimensional 
HB String Theory on AdSs x and four-dimensional J\f = A Super-Yang-Mills (SYM) with 
gauge group SU{N). This AdS 5 /CFT 4 duality is the most concrete example of the holographic 
principle and provides a non-perturbative definition of string theory. 

States of the field theory at large N and large t’Hooft coupling correspond to solutions of 
classical gravity in the bulk. In particular, bulk black holes describe thermal states on the field 
theory, with the field theory temperature T identified with the Hawking temperature of the AdS 
black hole. We thus expect black holes in AdSsxS^ to play an important role in understanding 
the phase diagram of A/" = 4 SYM. 

According to the correspondence, the background spacetime for the field theory is specified 
by the four-dimensional boundary of AdSs. Since J\f = 4 SYM is a conformal field theory, it 
does not exhibit phase transitions on a scale-invariant background like Minkowski space 
Instead, a different background spacetime can be chosen which allows for a more interesting 
phase structure. The phase structure of such solutions is a well-studied topic (see, e.g. the 
review and references therein). However, much of this study neglects effects set by the 
curvature scale of the S^. 

In this manuscript, we construct new thermal phases where the plays an important role. 
One of the most well-studied backgrounds for the field theory is the Einstein static universe 
Mt X S^. We will therefore be concerned with solutions that are asymptotically global AdSs x S^. 
These solutions must satisfy the type HB SUGRA equations of motion. With only the metric 
g and Ramond-Ramond 5-form F(b) = d(7(4) turned on, these equations are given by: 

Gmn = Rmn - ^FmpqrsFm^^^^ = 0, = 0, 


-^"( 5 ) = *-^( 5 ) ) ( 1 - 1 ) 









where we have imposed self-duality on the 5-form. Perhaps the most well-known solution to 
these equations is AdSs x which in global coordinates is described bjQ 
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where ej^iypardy^ A dy^ A dy^ A dy^ A dy^ and Cabcde^x^ A dx^ A dx^ A dx^ A dx^ are the volume 
forms of the AdS 5 ^.nd base spaces, respectively. L is the AdSs length scale and the 
radius. These are required to be the same by the equations of motion. 

We are interested in black hole solutions that are asymptotically global AdSsxS^. There is, 
of course, the Schwarzschild BH family (hereafter abbreviated as AdSs-SchwxS^ BH or simply 
the AdSs-Schw BH) with horizon topology S^xS^. This solution can be written as 


ds fdt -\~ d~ T dfl 3 -\~ L dfl^ , FpupcFT ^pupar-) ^abcde ^abcde i (1*^) 

2 ^ 2/^2 \ 

where f ='^ + J 2 - ■ 


When the horizon radius vanishes, this solution reduces to global AdSsxS^. 

Let us review the phase diagram of this family of solutions in a given thermodynamic 
ensemble. There are two ensembles that differ significantly: the canonical ensemble (at fixed 
temperature) and the microcanonical ensemble (at fixed energy). 

In the canonical ensemble, we must consider Euclidean solutions at fixed temperature T, that 
asymptote to S^xS^xS^, with the radius of the identified as the inverse temperature T~^. 
At a given temperature, there are up to three solutions in this family that compete: small black 
holes, large black holes, and thermal AdS (which corresponds to a thermal gas of gravitons). 
Computing the free energy from the Euclidean action determines which solution is preferred. 
At high temperatures, the large black hole phase dominates, while at low temperatures the gas 
of AdS gravitons is the dominant phase. Small black holes are never preferred. This indicates 
the existence of a critical temperature where a first order phase transition - the Hawking-Page 
transition - occurs [^[^. In the dual CFT 4 , this is interpreted as a confinement/deconfinement 
transition [^. 

In the microcanonical ensemble one fixes the energy and computes the entropy to find the 
dominant phase. In the AdSs-SchwxS^ family, there is only one solution at a given energy, so 
its phase diagram is trivial. 

However, the AdSs-SchwxS^ family assumes a large amount of symmetry. There may be 
other solutions with AdSsxS^ asymptotics which do not globally respect the SO{Q) symmetry 
of the SA 

Why would one expect such solutions to exist? As first observed by [5]-[^, small AdSs- 
SchwxS^ BHs are unstable to a Gregory-Lafiamme (GL) instability when the ratio be¬ 
tween the horizon radius and the radius, r^/L, is smaller than a critical value. Indeed, 
a topologically S^xS^ black hole with a large and small resembles a black brane, and 
systems with such a large hierarchy of scales are generically unstable to GL-type instabilities. 

In analogous systems, such instabilities typically contain a zero mode which indicates the 
existence of a new family of black holes. This new family often connects to other black holes with 
different topologies through conical transitions. For example, in the phase diagram of Kaluza 


Klein black holes that are asymptotically 


10-24 , there are S^xS^ black strings that 


^Capital indices M, N,... are d = 10 indices, Greek indices ... are AdSs indices, and small case Latin indices 
a,b,... are indices on the S^. 
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are GL unstable when the is small relative to the S^. The associated zero mode connects 
these solutions to non-uniform strings which themselves connect to topologically spherical 
black holes. 

As another example, asymptotically flat Myers-Perry BHs in d > 6 dimensions with a single 


non-trivial (but large) angular momentum also suffer from a GL-type instability 25 -^ (known 


in this setup as the ultraspinning instability). The zero mode in this case connects the Myers- 
Perry BHs to lumpy (a.k.a. bumpy or rippled) rotating BHs with horizons. Some of these 
lumpy solutions then connect to black rings with horizon topology x 26,32,33 . There 

are also zero modes corresponding to higher harmonics on the which yield other lumpy 


solutions that connect to multi-horizon solutions like black Saturns or di-rings |34| . 

In the present case, zero modes in AdSs-SchwxS^ BHs would lead to new black holes with 
the same S^xS^ topology, but break the SO{Q) rotation symmetry of the S^. We call these 
new solutions “lumpy” black holes. These lumpy black holes are expected to connect to other 
(possibly multi-horizon) solutions with horizon topology S^, or other products of spheres like 
S4 xS^. From a dimensional reduction, the is interpreted in the dual field theory as a number 
of scalar operators that respect an SO(6) symmetry. The GL instability is therefore associated 
with spontaneous symmetry breaking, where these operators develop non-trivial VEVs. In this 
case, we will argue that this phase transition is first order. 

Thus, there are many different solutions that compete in a given thermodynamic ensemble. 
Since these new solutions appear only for black holes that are sufficiently small compared to the 
S^, they are not expected to dominate the canonical ensemble, which favours large black holes at 
similar temperatures. However, in the microcanonical ensemble at small energies, we expect the 
solution with the most entropy to be a black hole with spherical topology S^. More specifically, 
the most entropic spherical black hole is expected to be the most symmetric, preserving the full 
50(4) symmetry of global AdSs well as the largest subgroup of 50(6), which is 50(5). We 
therefore focus on solutions which preserve the of AdS 5 ^ round within the S^. 

Let us write the metric of suggestively as 


^ 1-^2 


-h (1 — . 


(1.4) 


These are just the usual spherical coordinates with a redefined polar angle x — cos 6. Now 
we must isolate the perturbative modes that preserve the round S^. Any smooth perturbation 
about AdSs-SchwxS^ can be decomposed into a sum of perturbations of scalar, vector and 
tensor types. These types are defined according to how they transform under diffeomorphisms 
of the S^. The modes we are searching for appear in the scalar sector. In the coordinates of 
0 and ( |1.4[ ), these take the simple form Sqmn = where Y^^x) is a scalar 

harmonic on the and £ is its quantum number. If we focus on harmonics that only depend 
on the polar angle x (and hence preserve an S^), £ labels the number of nodes along x. 

The zero mode of AdSs-SchwxS^ that appears with the largest horizon radius is £ = 1. The 
linearised field equations reduce to a single ODE in the radial coordinate r and can be solved 
for any positive integer value of £. For £ = 1, we find that the zero mode appears at the critical 
horizon radius 

- 0.4402373 L, (1.5) 

as first found by Hubeny and Rangamani in j^; AdSs-SchwxS^ with are unstable. 

At this zero mode, we expect a family of lumpy black holes to emerge. By drawing an analogy 
with similar situations, we conjecture that this family leads to a conical merger with a family 
of black holes with topology, which can be thought of as being localised on the S^; see Fig. 

Eb. 
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Figure 1: A pictorial representation of the for some black hole solutions that asymptote to 
global AdSsxS^ (the last three are conjectured solutions). The first one is the familiar 
AdSs-SchwxS^ that is smeared over the S^. The second is a £ = 1 BH localized on 
the north pole of the S^. The third is a ^ = 2 localised BH solution and the fourth is a 
^ — 2 black belt. In this manuscript, we construct the £=1,2 “lumpy” BH solutions, 
with horizon topology S^xS^ that we conjecture to be connect to these localised BH 
solutions in a phase diagram. 


Of course, the £ = 1 zero mode is only the first mode that appears. There are an infinite 
number of such modes, with higher modes appearing at ever smaller horizon radii. For example, 
the £ = 2 mode appears at 

-0*3238898 L. ( 1 . 6 ) 

We note, however, that there are important differences between the even £ modes and odd £ 
modes. If 8g is a linear perturbation, then —8g is also a linear perturbation. In the odd £ 
modes, these can be mapped to each other via a Z 2 symmetry of the S^, and so are equivalent. 
For example, among the £ = 1 solutions, the choice of sign merely selects whether the localised 
black hole will develop on the north or south pole of the S^. In the even £ modes, however, 
these perturbations map to themselves under this symmetry. The 8g and —8g perturbations 
are not equivalent, which means we have two branches of solutions emanating from the even £ 
zero modes. In the £ = 2 modes, we expect one branch to lead to two disconnected black 
holes localised on the poles of the S^; see Fig. [l}c. We expect the other branch to lead to an 
s^xS^ black hole (the s^ being a smaller sphere than the S^); see Fig. [^d. The larger wraps 
around (coincides with) the equator of the S^, so we call these conjectured solutions “black 
belts” (the s^ gives the transverse directions of the belt). Higher £ modes lead to various other 
multi-horizon solutions with some combination of holes and s^xS^ belts. 

In this paper, we construct these lumpy black holes connected to the £ = 1 and £ = 2 zero 
modes and study their thermodynamic properties. We detail our numerical construction in 
section 2 , and compute the phase diagram in section 3. In appendix A, we give the technical 
details of Kaluza-Klein holography necessary to interpret our results on the CFT 4 (see also 
39,41.^)- Numerical checks are in appendix B. 
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2 Numerical Construction 


2.1 Lumpy solutions with ^ = 1 


Here, we give the numerical details of our construction of the £ = 1 lumpy black holes (BHs). 
We use the DeTurck method, which proceeds as follows. First, we choose a reference metric g 
that satisfies our desired boundary conditions (i.e. contains a regular horizon, has the correct 
asymptotics, and has the desired symmetry axes). Then, rather than solve the equations O), 
we instead solve the similar equations 


gS = Rmn - ^FmpqrsFm^"^^^ - ^(M^N) = 0, = 0, 


-^(5) = '^Fi 


(5) , 

( 2 . 1 ) 

where — g^^l^pgig) — ^pqig)] and r(gf) is the Levi-Civita connection associated with the 
reference metric g. 

The benefit of this method is that unlike the equations (0, the above equations 0 form 
a set of elliptic PDFs [^. For solutions of 0 to also be solutions of 0 , we must have 
— 0. In some cases (such as vacuum-Einstein), there is a proof that all solutions to B 
also have = 0. In our case, we do not have such a proof, so we must verify that = 0 
after solving the equations. The local uniqueness property of elliptic equations guarantees that 
solutions with ^ 0 cannot be arbitrarily close to those with = 0. 

Now we must find a suitable reference metric g. The most obvious choice is to use AdSs- 
SchwxS^, but rather than use the coordinates in (1.3) and (1.4), we instead redefine 

r+ 


r = 


1 


X 


= x'^ 2 — 


( 2 . 2 ) 


and set = Vp/L. With these new coordinates, the solution (1.2) can be written 

4y+ dy^ 


ds^ = 


„2'i2 


^4) = 


(1 -y 

+ L 2 

(2 - y 2 ) 

(1 - y2)4 


(2 - y^) G{y) dt^ + 
4da;^ 


(2 - y2) G{y) 
+ (l - a;^)^dfl|) , 


+ y+ df 23 


L4 


H (y) dt A 45(3) - ^ 45(4) , 


(2.3) 


where 


G(y) = (1 - y2)' + ylHiy) , H{y) = 2 - 2y^ + y\ (2.4) 

We choose this to be our reference metric g. Here, x G [—1,1], with x — —1 and x = 1 
corresponding to the north and south poles of the S^, and y G [0,1] with ^ = 0 corresponding 
to the horizon and ^ = 1 to the boundary of the AdSs. 

Using this reference metric, we write down a general ansatz which is static and preserves 
the symmetries of the and S^: 


d5^ = 


( 1 - 2 /^)^ 

-y^ (2 - y^) G{y) Qi dt^ + 

4dx2 


4y4 




Q4 


2 — x- 


(2 - y2) G(y) 

2 +Q6 {1-x‘^fdnl 


Q 2 [dy + (1 - y‘^)^Q3 dx] ^ + y+Qs dQj 


Ci4) = 


Ly^ y_(2— y_^ ^ -^w d 5 ( 4 ) 

(l-y 2 )^ 


(2.5) 
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where {Qi^W} (/ = 1,...,7) are functions of x^y which will be determined numerically. If 
we set Qi = (^2 = Q 4 = Qs = Qe = Qt = = 1 si^nd Qs = 0, we recover the AdSs- 


SchwxS^ solution as written in (2.3). On a solution different from AdSs-SchwxS^, our ansatz 


would preserve the full 50(4) symmetry of the and an 50(5) symmetry of the S^, allowing 
deformations along the polar direction x. Also, note that L drops out of the equations of motion, 
so we are left with one parameter and 8 functions {Q/, W}. 

At this point, let us count the number of equations we have. The tt, xx, x^, yy^ 
and 124124 components of the Einstein-DeTurck equation give 6 equations, the tl 23 and I 24 
components of the five-form equation give 2 equations, and the xl 24 and ^124 components of the 
self-duality condition give us an additional 2. With 10 equations and 8 functions, this seems 
like too many equations. Furthermore, the self-duality equations do not yield second-order 
differential equations, so our equations are not manifestly elliptic. 

It turns out that these issues do not pose a problem. One can use the two independent 
components of the self-duality constraint to algebraically solve for dxW and dyW. After sub¬ 
stituting these derivatives into the five-form ^( 5 ), one finds that F( 5 ) is now independent of 
VE, i.e. just a function of and their first derivatives. This means that differentiating 

to obtain the second derivatives of W and substituting those into the remaining equations of 
motion eliminates W completely. The five-form equations of motion reduce to a single equation 
and, together with the Einstein-DeTurck equations, we are left with 7 elliptic equations for the 
remaining 7 functions Q/. Since the metric and are now independent of VE, there is also 
no need to compute W to extract physical quantities. 

Now let us discuss boundary conditions (BCs). The BCs at the horizon y = 0 and the 
poles X = ±1 are determined by regularity. The conformal boundary ^ = 1 is determined by 
demanding that the solution is asymptotically AdSsxS^. More specifically, we would like the 
various operators in the dual field theory to be unsourced. Determining the correct powers of 
I — y that accomplishes this is a subtle issue which we defer to section A.5 of the appendix. 
Here, we simply give our boundary conditions in full: 


BCs at the poles (x = ±1) : < 


dxQi{y, ± 1 ) = 0 , 

Q3(2/,±1) = 0, 

Q4(?/, ±1) = Q6(2/,±1) , 


/ = 1,2,5,6,7; 
7 = 3; 

7 = 4. 


BCs at horizon (^ = 0) : 


Qi(0,x) = Q 2 ( 0 ,x) 
5yQ/(0,x) =0, 
( 53 ( 0 , x) =0, 


7 = 2,4, 5,6,7; 
7 = 3. 


( 2 . 6 ) 


(2.7) 


BCs at conformal boundary (^ = 1) : 


( 2 . 8 ) 


Qi(l,x) = l, 7 = 1,2,4,5,6,7; 

^ Q3(1,x) = 0, 7 = 3. 

Now we are in a position to solve the PDF system (|2.1|) subject to the boundary conditions 


(2.6)-(2.8). To solve the equations, we use a standard Newton-Raphson iteration algorithm 
based on pseudo-spectral collocation on a Chebyshev grid. We obtain a seed by perturbing 
the AdSs-SchwxS^ solution near the zero mode. Our solutions are parametrised by We 
will present our numerical results in section and discuss numerical checks in section [B] of the 
appendix. 


2.2 Lumpy solutions with 7 = 2 

It turns out that black hole solutions branching from modes with even values of 7 are easier to 
construct than those with odd values of 7. In order to see this we note the following two Z 2 
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symmetries: (1) modes with even t preserve the symmetry x —x of the line element (1.4); 


(2) they also preserve the symmetry r —r of the line element (1.3). While the first symmetry 
is easy to understancQ the second requires some explanation. 

It turns out that generic infinitesimal perturbations about AdS 5 xS^ can always be reduced 
to the study of a set of decoupled ODEs [^. Solutions to these ODEs can then be used to 
reconstruct metric perturbations, which in turn allows us to read off the decay of the several 
metric functions as we approach the boundary. Indeed, these decays are related to the mass 
of each perturbation from the AdSs perspective. It turns out that for even values of all 
perturbations lead to masses that cause the symmetry r ^ —r to be preserved. 

We can thus take advantage of these symmetries to design a better line element for a reference 
metric. We will still choose as a reference metric that AdSs-SchwxS^, but rather than use the 


coordinates in (2.2), we instead redefine 


r = 


Ly+ 


= x\/2 — 


(2.9) 


With these new coordinates, the solution (1.2) can be written 


ds^ = 


^4) = 


1 - y2 


-y^G{y) dt^ + 


y+ 


(1 -y2)G(y) 


+ y\ dfl 3 


+ L^ 


V2 (1- 


^ ^ ^ ^*^(3) - ^ d5(4) , 




4da:^ 

2-x^ 


+ (l — x^)^ dCll 


( 2 . 10 ) 


where 


G{y) = l-y^ + ylH{y), 


H{y) = 2-y\ 


( 2 . 11 ) 


We choose this to be our reference metric g. Here, x G [0,1], with x = 0 and x — 1 corresponding 
respectively to the Z 2 axis of symmetry and north pole of the S^, and y G [0,1] with ^ = 0 
corresponding to the horizon and ?/ = 1 to the boundary of the AdSs. 

Using this metric, we write down a general ansatz which is static and preserves the above 
Z 2 symmetries and also those of the and S^: 


d5^ = 




G4) = 


1 - y2 ^ 

+L2 

lSX 


~y^G{y) Qi dt^ + 


yl 


(l-y2)G(y) 


Q 2 [dy + (1 - y^)Q 3 dx] ^ + y+Qs df2| 


Zi X 


(I- 


H{y) Qt dt A d5(3) -^W 

(1 - y^) V2 


d5, 


(4) 


( 2 . 12 ) 


where {Q/, VE} (/ = 1,..., 7) are functions of x, y which will be determined numerically. The 
construction of these solutions and respective boundary conditions parallel those the previous 
section, and as such will not be presented here. The only change occurs at the new boundary 
X = 0, since there we must demand reflection symmetry. This is reduces to 


BCs at the reflection plane (x = 0) : 


dxQi{y,d) = 0 , 

<33(y,0) = 0, 


1 = 1 , 2 ,4, 5 ,6,7; 

1 = 3. 


(2.13) 


^Recall that the parity of scalar spherical harmonics under the symmetry x —x is simply (—1)^ 
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3 Phase diagram of AdSsxS^ BHs. 

In this section we wish to compute the phase diagram of the lumpy BH families and the AdSs- 
Schw BHs. We first need the thermodynamic quantities of the lumpy BHs. The BH temperature 
is proportional to its surface gravity at the horizon, 2ttT = [—with K being 


the generator of the Killing horizon. It follows from 
of the lumpy BHs and AdSs-Schw is then 


r = 


1 + 1 
L 27ry+ 


, with K = dt, that the temperature 


(3.1) 


The entropy of a BH is proportional to its horizon area, S = Hj^/(4Gio). From the AdS/CFT 
dictionary, the 10-dimensional (Gio) and 5-dimensional (G 5 ) Newton’s constants are related to 
N by 

^ Gio 


TT 




G 5 — 


7V^L^ 


The entropy of the ^=1,2 lumpy BHs can then be written as 


S = L^N^^yl 


dx ^^ \/Q 4 { 0 ,x)Q 5 { 0 ,x)^Qe{ 0 ,x)^, 


(3.2) 


(3.3) 


while the AdSs-Schw BH entropy is still given by this expression with Qd = Qs = Qe = 1 he., 

5schw = F3iV27ry3. 

Computing the energy of the lumpy BH is non-trivial. It can be read from the asymptotic 
expansion of the fields at the holographic boundary using the formalism of Kaluza-Klein (KK) 


holography and holographic renormalisation (see also [35}|39| 41,4^). A detailed discussion 
of this formalism applied to this system is given in Appendix]^ It culminates with the expression 
for the energy (A.54), that we reproduce here (this is valid for the £ = 1,2 lumpy solutions). 


E = 


7V2 


3072L2 


576 + 2304 yl{l + yp - yA5 /Ja + 30 /?! + 12 (16 So + ^ 4 ) 


(3.4) 


In this expression, {f32^So^64} are undetermined coefficients that appear in a Taylor expansion 


of the fields about the holographic boundary (see (A.7)-(A.IO)) after imposing appropriate 
Dirichlet boundary conditions (these correspond to having no sources in the dual CFT 4 ). These 
coefficients can be obtained by differentiating our numerical results. The energy of the AdSs- 
Schw BH is given by 

^schw = [(3/4)y2 (y 2 + 1 ) + 3/16] ^ 2/^2 ^ ^3 5 ] 

and is recovered when /32 = = ^^4 = 0. The AdSs-Schw BH energy contains a contribution 

from the AdSs background, E^dS^ ~ which is the well known Casimir energy of the dual 

A/' = 4 SYM on M X S^. 

Extracting the constants {/ 325 ^ 05 ^ 4 } requires accurately evaluating four or two derivatives, 
for odd or even respectively. Rather than evaluate up to four derivatives, a numerically 
simpler way is to integrate the first law dE — TdS. Some of our data has been extracted using 
this simpler method. Where we can accurately do both, these energies agree (this comparison 
is in Figs. [^and[^. 

With these thermodynamic variables, we can compute the (Helmholtz) free energy using 


E = E-TS. 


( 3 . 6 ) 
















For the AdSs-Schw BH this is given by 


^ V 3 1 2 2 N 

^SAdS^ ^ I 16 ^ 4^+ 


(3.7) 


When we set the horizon radius y+ = 0 we get the free energy of AdSs, TAdSz ~ 

Now let us discuss the various thermodynamic ensembles. For our system, these are the 
microcanonical and canonical ensembles. In the microcanonical ensemble, the energy is held 
fixed and solutions with higher entropy are preferred. In the canonical ensemble, the tempera¬ 
ture is fixed and the solutions with the lower free energy are preferred. To obtain adimensional 
quantities, we scale by factors of Gio ^ 1/A^^ and L. The relevant phase diagrams are therefore 
S/{L^N‘^) vs el ?for the microcanonical ensemble and TL ?vs TL for the canonical 
ensemble. 

Before giving the phase diagram with the lumpy BHs, it is instructive to first discuss the 
AdSs-Schw phases. In the microcanonical ensemble, the AdSs-Schw BHs have increasing en¬ 
tropy with increasing energy. As mentioned earlier, the zero entropy solution corresponds to 
the Casimir energy EAdS^ ~ At a given energy, there is only one solution in this family, 

so the phase diagram is trivial. 


AdSs-Schw BHs are more complex in the canonical ensemble. From (3.1), one can see that 
AdSs-Schw BHs have a minimum temperature at 


y+| = ^ ^ 0.707107, = — ^ 0.450158. (3.8) 

y2 tt 

There are thus two AdSs-Schw BH solutions with any given temperature above Tcy. These can 
be distinguished by their size, so AdS 5 -Schw BHs with horizon radius yj^ < yA are called 
small BH’s, while those with yj^ > yA are large. The free energies of large AdSs-Schw BHs are 
always lower than that of the corresponding small AdSs-Schw BH with the same temperature. 
That is, large BHs are preferred over small BHs. In the EL?/N‘^ vs TL phase diagram, the 
large and small BHs meet at a cusp at Tcy. There is, however, a third phase which is thermal 
AdSs xS^. This is just the Euclidean solution of AdS 5 xS^ with an arbitrary period chosen for 
the Euclidean time circle (i.e. at any temperature). Below a temperature corresponding to 

y+1 = 1 , ThpL 0.477465, (3.9) 

thermal AdSsxS^ has lower free energy than both large and small AdSs-Schw BHs. At tem¬ 
peratures above Thp^ large AdSs-Schw BHs are preferred. This is a first-order phase transition 
known as the Hawking-Page (HP) phase transition [^. In the dual CFT, this is interpreted as 
a confinement/deconfinement transition . 

This phase diagram will become more complex if we further allow for solutions that break 
the symmetries of the S^. As discussed earlier and in [^, the largest wavelength GL zero mode 
{i = 1) that preserves an SO{5) symmetry of the appears at 

- 0.440237, TqlL - 0.501653, (3.10) 

which is a horizon radius corresponding to a small BH. A branch of lumpy BHs emerges from 
this zero mode. 

Now let us consider these lumpy BHs in the microcanonical ensemble whose phase diagram 
is in Fig. Because the entropy between these lumpy solutions and AdSs-Schw are close, we 
have instead plotted the entropy difference between these solutions at the same energy. The 
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E/N^ 

Figure 2: Phase diagram in the microcanonical ensemble with the i = 1 lumpy BH family. We 
plot the entropy difference between each solution and that of the AdSs-SchwxS^ BHs 
as a function of energy. The red line represents the AdSs-SchwxS^ BH, while the 
blue dots describe the £ = 1 lumpy BH family. The green diamond marks the GL zero 
mode where the lumpy and AdSs-SchwxS^ families merge. We have computed the 
thermodynamic quantities both by integrating the first law (dark blue curve) and by 
reading the asymptotic energy (light blue curve; see text). The inset plot is a zoomed 
out version where we plot the entropy of the solutions as a function of their energy. 
Here, we fix L = 1. 


actual entropy is in the inset plot. There, we can see that these lumpy BHs always have lower 
entropy than AdSs-SchwxS^ BHs and therefore do not dominate the microcanonical ensemble. 

The phase diagram in the canonical ensemble with the i = 1 lumpy BHs is shown in Fig. 
The zero mode lies in the small BH branch, so we plot the difference in free energy between the 
lumpy BH and that of the small AdS 5 -SchwxS^ BH as a function of temperature. In the inset, 
we show a wider view of the actual free energy vs the temperature. We can see that the lumpy 
BHs have higher free energy than both large and small AdSs-Schw BHs as well as thermal AdS. 
They therefore never dominate this ensemble either. 

If the full 10-dimensional theory is dimensionally reduced to a theory on AdSs, fhe lumpy 
AdSs xS^ BHs are reinterpreted as five dimensional BHs with non-trivial scalar fields. From the 
perspective of the field theory dual to the AdSs, there are nonzero expectation values for scalar 
operators. This is spontaneous symmetry breaking. These expectations values can be computed 
using the tools of Kaluza-Klein holography , whose discussion we defer to Appendix (see 
also [35}j4^ ). We find that some scalar operators develop non-trivial VEVs. The expectation 
values of two of these scalar operators are plotted in Fig. |^as a function of the temperature. 
As expected, we find that these VEVs vanish at the zero mode (and for the AdSs-SchwxS^ BH 
family) and then their amplitude grows monotonically along the lumpy BH family as we move 
away from the zero mode. 

Now let us move on to the t — 2 lumpy solutions. As we have mentioned in the introduction, 
there should be two branches of black holes that emanate from the £ = 2 zero mode. We 
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Figure 3: Phase diagram in the canonical ensemble with the £ = 1 lumpy BH family. We plot the 
difference in free energy between each solution and that of the small AdSs-SchwxS^ 
BHs as a function of temperature. The red line represents the small AdSs-SchwxS^ 
BHs. Again, the green diamond marks the zero mode where the £ = 1 lumpy and 
AdSs-SchwxS^ families merge. We compute the thermodynamic quantities both by 
integrating the first law (dark blue curve) and also by reading the asymptotic energy 
(light blue curve; see text). The inset shows a wider range of the phase diagram and 
plots the difference in free energy between each solution and thermal global AdSs x 
as a function of temperature. The black square marks the Hawking-Page transition, 
where large BHs dominate at higher temperatures and thermal AdSsxS^ dominates 
at lower temperatures. We again set L = 1. 


conjecture that one of these branches connects to two BHs, so we call this the “double black 
hole (BH)” branch. We expect that the other branch connects to topologically s'^xS^ BHs with 
the wrapping around the equator of the S^, so we call this the “black belt” branch. 

The phase diagram of the £ — 2 solutions in the microcanonical ensemble is displayed in 
Figj5| Near the zero mode, the black belt branch extends towards higher energy, but with 
lower entropy than AdSs-SchwxS^. On the other hand, the double BH branch extends towards 
lower energy, but with higher entropy than AdSs-SchwxS^. Some of these lumpy BHs are 
therefore favoured over AdSs-SchwxS^ (but might not be the dominant phase of the ensemble). 
Furthermore, the double BH branch contains a turning point at a cusp. Thus, there can be two 
lumpy BHs for a given energy. 

Now let us continue with the £ — 2 solutions in the canonical ensemble shown in Fig. 
Near the GL zero mode, the black belt branch extends towards lower temperature and higher 
free energy than the corresponding small AdSs-SchwxS^ BH. The double BH branch extends 
towards higher temperature and lower free energy than the small AdSs-Schwx BHs. So also 
in this ensemble, the double BH branch of the lumpy BHs is therefore favoured over small 
AdSs-Schw xS^ BHs, but still have higher free energy than the large BHs. 

The £ = 2 solutions also have non-vanishing scalar operators whose magnitude grows mono- 
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Figure 4: Expectation values (O 52 ) (Left) and (O^s) {Right) of the dual operators to the KK 
scalar fields and (here we set L = 1) for the £ = 1 lumpy BH. We find that 
near the merger we have the fitting: (O 52 ) ~ ^4 (1 — T/Tc)^ with A ^ 0.0147zb0.0001 
and a - 0.966 ± 0 . 001 ; and (Ocs) ~ (1 - T/Tc)^ with B - -0.1643 ± 0.0006 and 

P ~ 1.4739 ± 0.0006. 


tonically away from the merger. This is much like the behaviour displayed in Fig. for the 
£ = 1 case, so we do not present a separate plot. 

Putting the i = 1 and i = 2 plots together, we get Fig. [^for the microcanonical ensemble 
and Fig. for the canonical ensemble. We note that at the GL zero modes, the slope of the 
entropy and free energy of the lumpy solutions match that of the AdSs-SchwxS^ BH. This 
indicates a second-order phase transition and is consistent with the fact that these phases arise 
perturbatively. 

Now we attempt to analyse the approach of the lumpy solutions towards the conical mergers. 
Let us first discuss the £ = 1 case. In the left panel of Fig. we plot the Ricci scalar of the 
induced horizon geometry on each of the poles of the as a function of temperature. We see 
that the Ricci scalar is getting large at one pole and small at the other pole. In the right panel 
of Fig. 1^ we plot the radius of the at the horizon as a function of the polar variable x of 
the for four different temperatures. As the temperature decreases and we move away from 
the GL zero mode we find that this radius is getting small at the South pole (x = 1), consistent 
with the conjectured conical merger. As we mentioned earlier, we suspect topologically BHs 
on the other side of this conical merger. 

Now we proceed with the i = 2 case. In Fig. we plot the Ricci scalar of the induced 
horizon geometry at one of the poles of the and at the equator, both as a function of 
temperature. The curvature of the black belt branch is getting large at the poles. In Fig. [E 
we plot the radius of the on the horizon as a function of the polar angle x of the for 
four different temperatures. In the right panel, we see that close to the conjectured conical 
merger, the radius of the black belt branch gets very small on the poles {x — 1) of the S^. In 
the left panel, we see that the radius of the double BH branch is decreasing on the equator 
{x = 0 ) of the as we approach the conical merger, though we are still somewhat far from 
this conjectured merger. This also explains why the induced Ricci scalar is not yet appreciably 
large at the equator for this family (see Fig. 10). 
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Figure 5: Phase diagram in the microcanonical ensemble similar to Fig. but now with the 
£ = 2 lumpy solutions. The red line with AS{E) = 0 represents the AdSs-SchwxS^ 
BH. Black inverted triangles describe the £ = 2 double BH branch while the brown 
triangles represent the black belt branch. The magenta square marks the £ = 2 zero 
mode where these lumpy BHs and AdSs-SchwxS^ families merge. The inset plot is, 
again, a zoomed out plot. 



T 

Figure 6: Phase diagram in the canonical ensemble similar to Fig. but this time with the 
£ = 2 lumpy solutions. The red line with A£F{T) = 0 represents the AdSs-SchwxS^ 
BH. Black inverted triangles describe the £ = 2 double BH branch while the brown 
triangles represent the black belt branch. The magenta square marks the £ = 2 zero 
mode. Like in Fig. [Slthe inset plot gives a broader view of the phase diagram. 
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Figure 7: Phase diagram in the microcanonical ensemble that collects the information displayed 
both in Fig. [ 2 ] and Fig. The red line represents the AdSs-SchwxS^ BH; blue dots 
represent the £ — 1 lumpy BH family; black inverted triangles describe the 1 — 2 
lumpy double BH branch; and the brown triangles represent the black belt branch. 


CM 
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Figure 8: Phase diagram in the canonical ensemble that collects the information displayed both 
in Fig. [ 3 ] and Fig. The red line represents the AdSs-SchwxS^ BH; blue disks 
represent the 1—1 lumpy BH family; black inverted triangles describe the 1 — 2 
lumpy double BH branch; and the brown triangles represent the black belt branch. 
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Figure 9: Left Panel: Ricci scalar of the 1=1 lumpy BHs evaluated at the horizon and at the 
north (x = 1; upper curve) and south poles (x = — 1; lower curve) of the S^. Right 
Panel: Radius of the evaluated at the horizon as a function of the polar angle x 
of the for the 1=1 lumpy BH. The green diamonds are for the solution closer 
to the GL merger (T = 0.50167), while the blue dots describe the lumpy solution 
with the lowest temperature (T = 0.49444) we have reached. In between we have two 
other curves with intermediate temperatures, namely T = 0.50120 (empty squares) 
and T = 0.49898 (circles). These solutions appear to be approaching a localised black 
hole. 



T 

Figure 10: Ricci scalar of the induced horizon geometry for the i = 2 lumpy double BH branch 
(black inverted triangles) and i = 2 lumpy black belt branch (brown triangles). We 
show this quantity evaluated both at the pole {x = 1] filled triangles) and at the 
equator (x = 0; empty triangles) of the S^. 
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Figure 11: Radius of the evaluated at the horizon of the 1 — 2 lumpy solutions as a function 
of the polar variable x of the S^. Left panel: 1 — 2 lumpy double BH branch with 
the magenta squares being the solution closer to the GL merger (with T = 0.59448), 
and the black filled triangles being the solution with the highest temperature (T = 
0.60849). In between we also present the solutions with T = 0.59922 (circles) and 
T = 0.60552 (empty inverted triangles). Right panel: 1 — 2 black belt branch with 
the magenta squares being the solution closer to the GL merger (with T = 0.59448), 
and the brown filled triangles being the solution with the lowest temperature (T = 
0.57278). In between we also present the solutions with T = 0.59015 (circles) and 
T = 0.58283 (empty triangles). 
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4 Discussion and Prospects 


Let us now summarise our findings. As predicted by [^|^, the authors of found that small 
AdSs-Schw xS^ BHs suffer from a Gregory-Laflamme instability. This instability contains zero 
modes where new stationary solutions are expected to exist. We constructed the “lumpy” 
solutions corresponding to the £ — 1 and 1 — 2 modes that break the SO{Q) symmetry of the 
down to 50(5). Due to the symmetries of the linear perturbations, the ^ = 1 mode contains one 
branch while the 1 = 2 mode contains two branches. We called the 1 = 2 branches the “double 
black hole” branch and the “black belt” branch in anticipation of the ensuing (conjectured) 
conical mergers. In the microcanonical ensemble, only parts of the double BH branch of the 
1 = 2 mode is preferred over AdSs-SchwxS^ BHs, but we do not have any evidence that it 
actually dominates the ensemble. None of these phases are preferred in the canonical ensemble, 
which is likely dominated by large AdSs-Schwx BHs. 

We have good numerical evidence that these solutions approach conical mergers. We also 
have evidence that the i = 1 lumpy BHs transitions to a topologically BH that sits on one 
of the poles of the (see Fig. 0i>) . We also have evidence that the 1 = 2 double BH branch 

transitions to two BHs, localised on each of the poles of the (Fig. [^c), and the black belt 

branch transitions to an s^xS^ black hole with the smaller s^ wrapping around the larger 
equator of the (Fig. [^d). 

Let us now speculate on the complete phase diagram which we conjecture to be something 


resembling Fig. 12 Consider the microcanonical ensemble {Left Panel) and solutions with i = 1. 


Small localised BHs should look muck like d = 10 small asymptotically flat Schwarzschild BHs 
whose entropy scales as 5 ^ We can compare this with that of small AdSs-SchwxS^ BH, 

whose entropy scales as 5 ^ We therefore conclude that the entropy of a small localised 

BH is larger than the entropy of a small AdSs-SchwxS^ BH. We further expect these BHs to 
merge with the lumpy BHs with horizon topology 5^xS^ at some conical merger point (point 
B in Fig. 12). A simple curve that satisfies these two properties is the dashed blue curve CB 

[T^ (the turning point of this curve and similar curves in 


sketched in the Left Panel of Fig. 
the diagram must be a cusp to be consistent with the first law of thermodynamics). If this 
conjectured curve turns out to be correct, there will be a first order phase transition where the 
entropy of the localised BHs and the AdS 5 -SchwxS^ BHs are equal at the same energy. Such a 
phase transition can be interpreted in the dual field theory as spontaneous symmetry breaking, 
though here this would be a first-order transition rather than second order. 

Consider now the 1 = 2 solutions. Both the double BH and black belt should be unstable 
to the formation of a single localised BH. Indeed, the double BH is an unstable equilibrium 
configuration and a small perturbation should make the two BHs merge into a single localised 
BH. Similarly, if we slightly perturb the black belt along the polar S^ direction, it is reasonable 
to expect that it will also collapse into a localised BH at the pole. Another possible instability 
mechanism is the fragmentation of the black belt into one or an array of BHs along the equator 
of the S^, and these should again be unstable and merge into a single localised BH. For these 
reasons, we expect that the double BH and black belt should be less preferred phases than the 
£ = 1 localised BH, both in the microcanonical and canonical ensembles. The simpler scenario 


with these properties is described by the dashed black (CE) and brown (CF) curves in Fig. 12 
Consider now the canonical ensemble {Right Panel of Fig. 12). Again, small localised 


BHs resemble 10-dimensional asymptotically flat Schwarzschild BHs whose free energy scales 
as F ^ T~^. Small AdSs-SchwxS^ BHs have a free energy that scales as F ^ T~^. Thus, for 
high temperatures, the localised S^ BHs are likely to have a lower free energy than small AdSs- 
SchwxS^ BHs. Though, large AdSs-SchwxS^ BHs have a free energy that scales as F ^ T^, 
and they are perturbatively stable so they are likely to dominate the canonical ensemble. 
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Figure 12: Conjectured phases diagram of asymptotically AdSsxS^ static BHs in the micro- 
canonical {Left Panel) and canonical {Right Panel) ensembles {AF = F — FAdSPi 
L = 1). The red line is the AdSs-SchwxS^ BH, the blue line AB describes the i = 1 
lumpy AdSsxS^ BH family with horizon topology S^xS^. The green diamond A is 
the i = 1 zero mode. The blue dashed line BC represents the conjectured loealised 
AdSsxS^ BHs with horizon topology S^. We expect there to be some point B which 
is the conical merger between these solutions. The turning points must be cusps to 
be consistent with the first law. The black square is the Hawking-Page (HP) critical 
point. The magenta square D is the 1 — 2 zero mode with the double BH branch 
along the brown curve and the black belt branch along the black curve. There are 
conjectured points E and F that mark conjectured conical mergers. 


We should mention that, in the sketched phase diagrams of Fig. we are probably 
oversymplifying the structure of the solutions near the conical mergers. Indeed, it might well be 
the case that the lumpy and localised branches will spiral towards the conical merger, leading 

In the present 


to an infinite discrete non-uniqueness similar to the one found in [33,34,43 45 


case, we do not approach the conical mergers sufficiently enough to address this question. 

Of course, to fully complete these phase diagrams, the localised solutions need to be con¬ 
structed. We leave this to future work that is currently in progress [^. We note that it is not 
necessary to resort to numerics to contribute to our understanding of this phase diagram. In 
particular, small localised BHs and black belts should be well described by black branes, and 
are hence amenable to a matched asymptotic expansion or a blackfold approximation (similar 
to the analysis done in (T5}|l7[[^[^[47[[4^ for localised BHs on a 46 1)0 

We also note that we have studied but two modes in the entire spectrum of spherical har¬ 
monics on S^, and we have only focused on those preserving an 50(5) symmetry. The full phase 
diagram is thus incredibly rich. Though, since the localised BHs connected to the i = 1 modes 
would possess a full 50(4) x 50(5) symmetry, they are likely the most symmetric of the single 
localised BHs, and are thus likely to be the entropically dominant phase for small energies 
in the microcanonical ensemble. For £ > 1, we can have multi-BH configurations localised on 
the even with different sizes. 

We have also only focused on global AdS^ which corresponds to a field theory background 
on M X S^. Other backgrounds such as x 5^ 49 or BH backgrounds yield gravity 


solutions with physics near the AdS scale. It would be interesting to understand how breaking 
the symmetries of the will influence these geometries. 


^We note that although asymptotically flat BHs only have a blackfold description when one of the spheres is odd, 
the s^xS^ belts are not supported by angular momentum, but by the geometry of the S^. 
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Our choice of boundary conditions ensures that the CFT dual to our gravity solution is 
M — ^ SYM [^. However, there are other choices that are consistent with finite energy and the 
absence of ghosts. Most notably, we could have chosen boundary conditions that correspond 
on the CFT side to adding a relevant double trace deformation to A/" = 4 SYM [^. The effect 
of these boundary conditions on the phase diagram remains unclear. 

There are also a number of other known AdS/CFT dualities such as those arising from 
AdS 4 X and AdSr X S^. Where the asymptotics are global AdSq x S^, we expect similar 
behaviour to what we have found in AdS 5 x S^. 
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A Kaluza-Klein holography 


In the main text, we have omitted any details concerning Kaluza-Klein (KK) holography [40] 
which we will supply here in this appendix. There are three tasks that required the use of this 
KK formalism: computing the energy, computing the VEVs of the scalar fields, and determining 
the appropriate asymptotic boundary conditions that correspond to turning off sources on the 
field theory. 

The gauge invariant formalism of KK holography was developed by Skenderis and Taylor 
in 1^. Previous studies useful for this endeavour are [35}j^ , and KK holography is further 
discussed and applied in 41,^. We will review this formalism in some detail, following (40|. 


We will also need to extend the results of 40 to our system where odd harmonics are excited. 


Moreover, we use a different harmonic representation for the S^. We will try to be self-contained 
but refer the reader to [^ and 35,^ for a more thorough exposition. 


The aim of KK holography is to first dimensionally reduce solutions with AdS^ x asymp¬ 
totics to solutions on AdS^, then apply holographic renormalisation to compute field theory 


quantities on the boundary of AdS^ 


40 


In our case (asymptotically AdSsxS^ solutions), the 


dimensional reduction requires expanding any solution as a sum of harmonics of the S^. These 
harmonic modes are interpreted as fields in the reduced AdSs theory. The behaviour of these 
fields on the boundary of AdSs give VEVs of operators on the dual conformal field theory. 

The dimensional reduction obtains the effective d — b fields T from some d = 10 fields 
In general, the map between T and ^ is highly nonlinear. However, if we are only interested 
in computing VEVs in the dual field theory, we only need the field T up to some order in a 
Fefferman-Graham expansion off the AdSs boundary. That is, for any particular VEV, we can 
write T as some polynomial of ^ and its derivatives, truncating at a particular order jioj. 
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For example, the quadratic expansion for a field takes the form 

= + + 0 {[^ 

Im 


(A.l) 


for some constants Jkim stnd Lkim- If is dual to an operator of dimension /c, then we would 
require expanding off the boundary in a Feffermann-Graham expansion to 0{z^). The 
quadratic terms with I -\- m — k also contribute to such an expansion. Higher order terms 
contribute as well, but for a given fc, we can truncate this expansion [^. In our case, it suffices 
to stop at quadratic order. 

At this point, it would be useful to give a brief overview of the lengthy calculation to follow. 
We begin in section A.l by writing our lumpy BH solutions as deformations of global AdSsxS^, 
and expressing those deformations as a sum of harmonics We then carry out a 

Feffermann-Graham expansion of the coefficients of these harmonics ^ to the order needed to 
extract VEVs. 

To avoid gauge issues, we now need to rewrite these coefficients in a gauge invariant way, 
keeping up to quadratic terms in the number of fields (recall we only need up to second order 
to extract the VEVs of interest). In section ( |A.2 ), we write down the gauge invariant quantities 
at linear order in the number of fields (and call these collective quantities t/^), and also show 
that most of these gauge-invariant fields obey an effective Klein-Gordon equation for a massive 
scalar. We proceed with quadratic order in section A.3 , and call the resulting quantities The 
second-order fields obey an inhomogeneous Klein-Gordan equation for a massive scalar. The 
source term depends on the square of linear-order fields and their derivatives. 

We perform the KK reduction in section[A]4, There, we w ill ob tain the effective 5-dimensional 

we will obtain the effective 5- 


A.5 


field in terms of the 10-dimensional ^ and In section 
dimensional action, namely (A.45), which describes KK scalars subject to a certain potential 
and living in the 5-dimensional background G with a negative (5-dimensional) cosmological 
constant. Section A.5 ends with the Einstein equation (A.47)-( A.48[ ) that the reduced graviton 
G obeys. This equation has a non-trivial energy-momentum tensor. 

Section A.5 also takes the 5-dimensional gravitational {G) and scalar fields (4^) and applies 

We first introduce the Fefferman- 


the standard holographic renormalisation procedure 
Graham coordinate Z — Z{z) for the 5-dimensional solution, and then do the standard Fefferman- 
Graham expansion off the 5-dimensional AdS boundary Z = 0. We can then construct the 
associated holographic stress tensor (A.53) and VEVs (A.51) of the most relevant KK scalars. 
In particular, from the holographic stress tensor we can then read off the expression for the 
energy, which we use in the main text. In the process, we also explain our physical motivation 
for our choice of asymptotic boundary conditions. 

Summarising our notation, describes the coefficients of the 10-dimensional harmonic ex¬ 
pansion around global AdSsxS^; represents gauge invariant quantities at linear order in the 
number of fields; p; describes gauge invariant quantities at quadratic order; and finally 4^ de¬ 
scribes the reduced 5-dimensional KK field (i.e. capital letters denote 5-dimensional fields and 
lower-case letters always refer to 10-dimensional fields). 


A.l Lumpy AdSsxS^ BHs expanded in spherical harmonics of the 

Our lumpy BH solutions are asymptotically AdSsxS^. Before proceeding with a dimensional 
reduction to AdSs, we first need to expand these solutions in terms of harmonics on the S^. 
There are scalar, vector, and tensor harmonics, which are defined by their transformations on 
the S^. The details of the harmonic expansion differ for each type, but only the scalar harmonics 
are consistent with our preserved symmetries (a SO{5) subgroup of 50(6)), so we will only focus 
on these harmonics. 
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If the is written as 


dnl 


AdX^ 

2-X2 


+ {1-x^fdnl, 


the regular (axisymmetric) scalar spherical harmonics are given by 


Yi{X) = 


{£ + 2){i + 3) 

2i(^+i)V3 


(^-A^ + 4;|i(l + Xy2-X2)^ , 


(A.2) 


(A.3) 


and satisfy 


Us^Y(,{X) = K(,Yi{X), with A^ = -£(£ +4), ^ = 0,1,2,. 


(A.4) 


The quantum number ^ is a measure of the number of nodes along the polar direction X that 
was quantised by requiring regularity at the poles X = ±1 of the S^; we set the azimuthal 
quantum number m = 0 because these modes would further break the SO{b) symmetry. We 
have chosen a normalisation in (A.3) so that 


with 


/s^ 


4^) = 


O 5 


+ !)(£ +2) ’ 


fis = TT^. 


(A.5) 


Now let us expand our lumpy BH solutions in terms of these harmonics. First, we write the 
fields of the solution as a deformation of AdS 5 x 


gMN = 9 mn Y hiMN, ( A . 6 ) 

Fmnpqr = Fmnpqr + Imnpqr , 

where {g°,F°^y is global AdSsxS®. Here, h and / need not be small. The field fluctuations 
abound global AdSs x S® thus admit the harmonic expansion: 


h^yz,X) = J2hUy)YiiX), 


i 


h^a{z,X) = J2^y)DaYi{X), 


e 


h(^ab){z,X) = Y,¥{z)D^aDb)YyX)) , 

i 


K{z,X) = J2^yz)Ye{X), 

i 


(A.7) 


and 


f/j,iypar{^-) X) = J25D[^bi^^^yz)Ye{X), 

i 

fab/iupi^’) — 0 5 fabcpui^’) — 0 5 

fabcdp{z,X) = eabcdD,Y^{X) , 

I 


fa^.upa{z,X) = Y,bUpa{z)DaYyX) , 

(A.8) 

fabcde{z^X^ E byz)AieabcdeYi{X), 


where we use the symmetric traceless notation = ^{Aab + A})a) — \gabA^. It follows from 
the field equations that there is an algebraic relation between the coefficients Uypcrr h^. 
Therefore we do not discuss the fluctuations &5pcrr further. Note that the case £ = 1 is 
special because = 0 so is not defined. The case £ = 0 is also special since A^ = 0, 

DaY^ — 0 and D(^aDb)Yi — 0; therefore, 0^^^, and are not defined. The expansion of 

all other fields start at £ = 0. 

In anticipation of reading off field theory quantities from the boundary of the reduced AdSs, 
let us also perform a Feffermann-Graham expansion off the boundary. The “background” fields 
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that give AdSsxS^ are expanded as 


L2 


dso = 


dz^ 


F° — f 

^fii^par ~ ^pvpcrTi 


16 

d^abcde ~ ^abcde i 


l + :,z^ + TzZ^ + 0{z^)]dt^ + {l--z^ + -z^ + 0{z 
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(iflq 


+ L^^dnl 
(A.9) 


where we have stopped at 0{z'^) which contains terms necessary to compute the holographic 
stress tensor and the VEVs of scalar operators. Note that the factors of 2:^/2 are present 
because the background is global AdSs, planar AdSs, the factors of z^/lG come from the 
conformal anomaly of AdSs. We have chosen a conformal frame where the boundary geometry 
of the AdSs is M x 5^. 

Now let us proceed with the first few harmonic coefficients in (A.7) and (A.8). Performing 
a Taylor expansion of the lumpy BHs up to 0(^^), or {L‘^/z‘^)0{z'^) for fields on the AdSs base 
space, the equations of motion yield (henceforward, we set L = 1): 
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5^7 

36864 


TT = — 


A 
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/32 ( 642;2 + (l7y^ P 2 + 20) z^ 


TT = 


y+ 


732 ; , 


6^=1 = 0, 


2 

6^=2 = 


y+ 


10 1024 ' 


P 2 ( 160 + (31 yX P 2 + 60) z" 


° “^1/7 73^ , 


y+ 


Vfdi' 5 


yl (34/3| + 5/32 - 20 ^ 4 ) ; 


^ V 7 y| (34^1 + 5/32 - 2054) 2;^ : 


(I16/3i + 25/32 - 10054) 


(A.IO) 


where = diag{ —1, 0, rjij} with rjij being the line element of a unit radius on the AdSs base 
space. To shorten the presentation, we have introduced the auxiliary constants {C'w? C'o? C' 2 } 
such that {Cw, ^^ 0 ,(^ 2 } = {1,1,1} for /i = u = t, while {Cw, ^^ 0 , 02 } = { — 1/3, —1/3,1/9} for 
fi = u = x'^) in and These harmonic coefficients depend on the horizon radius 

and on four undetermined constants {/S 25 73 , ^ 4 } that appear in the Taylor expansion off-the 

boundary z = 0. In this expansion, we have already imposed appropriate asymptotic boundary 
conditions (BCs) that eliminate several extra undetermined constants that would appear in the 
expansion. We will defer our discussion of these BCs to a later section A.5 ; see in particular 
the BCs (A.50) and Fig. 13 therein. 
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For reference, the leading terms of the transformation that bring the lumpy BH from the 
{x^y} coordinates of our ansatz (2.5) into the FG coordinates {z^X} are 


y = l-y±z-y^z^-^y+{2A + yl [l2 +/32 (l - 6X^(2 - X^))]) + •• • , 

a; = X-^y2 ^2X(2-3X2 + xV 

+ ^y2-X2 (1 - X2) yl [371 + 273 (X^ - 2X2 _ 1)] ^3 ^ . (A.ll) 


A.2 Gauge invariance and field equations at linear order 


Not all of the coefficients in ( |A.7[ ), ( |A.8| ) are independent. 

X 


Under a gauge transformation 


X + the fluctuations transform SiS h ^ Sh and / ^ 5/, where to linear order in the 
number of fields, these are given by 


ShMN = {Dm^n + Dn^m) + {DMi^hpN + Djsii^hMP + DpHmn)] (A. 12) 
S/mNPQR = ^D[m^^FnPQR]S + InPQRIS + ^^DsfMNPQR)- 

The gauge parameter itself has a scalar harmonic expansion, 

Uz,X) = , Uz,X) = J2^'H^)DaYe,{X). (A.13) 

4 4 


To avoid further gauge issues, we would like to use gauge-invariant quantities, which in this 
section we do to linear order in the number of fields. At linear order in the fluctuation, only the 
leading terms in the gauge transformation ( |A.12 ) contribute and the coefficients in the harmonic 
expansion (A.7) and (A.8) transform as (as justified above we do not need to discuss bjj^jjpcr) 


= + for ^>0; + for ^>1; (A. 14) 

for I > 2 - Stt^ = 2A^et for £ > 0; S¥ = for £ > 1. 


As we have mentioned earlier, some field coefficients are not defined for £ = 0 and/or I — 1. 

It follows that for ^ > 2 one can define three gauge invariant quantities: 

• For i >2 : — A^0^ , 

_ pi _ p 

¥ = (A. 15) 

with auxiliary field = =7 

For ^ = 1, since is not defined, the system is effectively described by one less gauge invariant 
quantity. One can define the gauge invariant quantities: 

• For TT^ = 0 , 

4^ = 4^ “ , 

6‘=6‘-^7r‘; (A.16) 

with auxiliary field ^ 6Bj^ — 

Z iV]^ 
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The case ^ = 0 is special since (j)^^ and are not defined. One also has ^ = 0 since 
= 0. It follows that is itself already gauge invariant and is not defined. This 
leaves the gauge invariant quantities: 


• For £ = 0 : tt^ = ^ 



iP not defined. 


(A.17) 


Note that is just a deformation of the background metric and the gauge invariant 

combination was chosen because it obeys the linearised field equations. 

With these gauge invariant quantities at linear order, one can introduce the gauge invariant 
combinations 




20(^ + 2) L 


TT^ - 10(£ + 4)6^ 


‘'= 20 (^ 

(A.18) 


with inverse relations b^ = —s^ + i^ and = lOks^ +10{£ + 4:)i^. To leading order in the number 

of fields, these satisfy the equations of motion 

as^ ^£{£-4)s\ for ^>1; □ = (^ + 4)(£ + 8) for £ > 0, (A.19) 

where □ is the D’Alembertian in AdSs. That is, the scalar field has mass = £{£ — 4) (in 
AdSs radius units L = 1) with the conformal dimension of the dual operators being {A+, A_} = 
{£,4 — -£}, while the scalar field F has mass = £{£ + 4)(£ + 8) which corresponds to the 
conformal dimensions {A+, A_} = {£ + 8, —^ — 4}. 

The massive KK gravitons also couple to the scalar harmonics (actually, the instability of 
the AdSs-Schw BH occurs in the KK graviton sector with £ > 1) and are described by the 
transverse and traceless fields 

iU = hi.) - + d*" - 12»') •'><>■ (A.20) 


which obey the equation 


□ 


<^U = [^(^ + 4)-2]<^U’ ^>0- 


(A.21) 


Thus, the KK gravitons (spin 2) have mass = £(^ + 4) and the conformal dimensions of the 
dual operators are {A+, A_} = + 4, —£}. For £ — 0^ the combination in (A.17) obeys 


the 5-dimensional linearized Einstein equations (the shift by rr^ can be traced back to the Weyl 
transformation required to write the 5-dimensional action in the Einstein frame). 


A.3 Gauge invariance and field equations at quadratic order 


At quadratic order in the fluctuation, all the terms in the gauge transformation (A. 12) con¬ 
tribute. For our purposes, it is not necessary to discuss the KK gravitons with £ > 0, so 
for compactness we omit these fields from our discussion. The interested reader can find the 
analysis of these modes in the original papers. 


24 









To discuss the gauge transformations, it is necessary to first project the fields into the 
spherical harmonic basis. This projection introduces the following integrals of the spherical 
harmonics {z(£) was introduced in the normalisation (A.5)): 

D(^aDb)Yi^D(^aDb)Yi^ = , with q{t} = !){£ + Y){^ + 5), 

ai23= [ Yi^Ye.Ye,, 6123=/ Yi^DaYe,D‘^Ye,, 

Js^ Js^ 

0123=/ D^^D^^Yi^DaYi.DbYi,, dus = [ Ye,D^^D^^Ye,DaDbYe, , 

Js^ Js^ 

6123 = {2DaD^Yi,D^,Db)Yt, + D^Yt,D,D^^,Db)Yt,) . (A.22) 

Up to second order, the coefficients in the harmonic expansion (A.7) and (A.8) transform 
as (recall that 

• For £i>2\ 

^ 5] [2y^e^^i^i23 + 


5^1 = + 




q{£i) z{£i 


^2/3>l 




V 1/ 0 ^ 


(A.23) 


^2/3>l 


• For £i — 1 \ 

5 ^^ — 2Ai(^^, 

5^^ not defined, 

5h^ (A.24) 


• For £ = {)'. 

= E ^ +2e,B^^)Ai 

i>i ^ ^ ^ 

5^^ not defined, 

dip not defined, (A.25) 



For notational convenience, in these expressions and all the expressions of this section, whenever 
we have with £ = 1 we mean This means, e.g. that in ( |a:231 ) 

is a short-hand notation for = E£2=2/3>1 ^7^^"^123 + E£2>2/3>1 


Similarly, one uses the short-hand notation E^>i = l^C^^9(l)+E£>2 

in (TA^. 
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Using these transformations one can check that the gauge invariant quantities to quadratic 
order areE] 


• For ii>2\ 




2^(4) 


E [(gA&a,23 + 6123 - 1 ^* 12 ) + (<il23 - ^6323 




+ 2^123 - 5g(£^) ^213 ] ) 


+2B 


f,2 I £>^*71^3 




ai23 + B^^{bi23 — ^^^Ci23) 




V8A,,"“' ' 20<i(<i) 




1 




1 

A^i 


Sf„ci23 + —D%H: 


'213 


^ C^213 

(A.26) 


• For ii = 1 : 

TT^l _ 

=V\ (A.27) 


• For 2^1 = 0 : 


^0 

TT = TT 


i>l 




1 

4 ' 


not defined, 


(A.28) 


where the linear order quantities tt^ and are defined in (A.15)-( A.16[ ). 

Although we do not discuss the details of massive KK gravitons (i.e. with ^ > 1; see 
Eol), the properties of the massless KK graviton will be fundamental for our later analysis. 
Under a gauge transformation this field transforms as 


sh% = [D^e+D.e^] + 




jiu 


(A.29) 


z{t) 


+ J2W - Uehl, + 2D^^eBi) 

i>l ^ 


and thus the corresponding gauge invariant massless graviton is given by is the AdSs 
metric) 


hi = hi + Inl;, - E ^ (-^A,(//^^. + \d,cI>^D^cI>^) 


£>1 


(A.30) 


+D^BJhl, + + BJDX.. + D,B^D,B^, + B^B^u - ) , 


'‘To get (A.26), symmetrise (A.23) and promote —)■ 4> —>■ Bf^} as follows from (A.14) and (A.15). 
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With these gauge invariant quantities at quadratic order one can introduce the gauge in¬ 
variant combinations 


/ = 


t^ = 


20(1 + 2) I 


TT^ - 10 {£ + 4)6^ 


for £ > 1 , 


1 


20(£ + 2 ) 

which obey the equations of motion 


(n^ + 10 £b^'j, for £> 0 , 


(AM) 


(□ - ml,) / = 
(□ - ml) = 


1 


(£ + A)(£ + 5)Qi + Q 2 + (£ + 4:)(DnQ^ + Q\) 


2 (£ + 2 ) 

1 

2(^ + 2) L' 


£(£-l)Q{ + Qi-£(D^Q^^ + Qi) 


(A.32) 


where □ is again the D’Alembertian in AdSs. The quantities {Qi^ ^ Q 4 } in the RHS of 

these two equations are the same as those defined in equations (3.27)-(3.34) of [^. On the 
other hand, the masses of these scalar fields are the same as that of the linear fields and 


mil = i{£ - 4), = {i + 4)(£ + 8). 


(A.33) 


Letting = { 5 ^, t^} and and using the equations of motion for CO^ 

and its derivatives one finds that the equations of motion (A.32) can be put in the form 


E [ 

A/3>1 

(A.34) 

for some coefficients 

A.4 5-dimensional KK description of the lumpy BHs 

Now we are in position to obtain the KK map. One can remove the derivative terms on the 


RHS of (A.34) by a field redefinition -ipf —)■ 4/^. The inverse of this relation gives the the reduced 
d = 5 field 4^^ = {S^,T^} in terms of the d = 10 field ip£ = {s^,t^} and is given by 


= w(ipei] 


- E(* 




(A.35) 


with the normalisation factor w{^pi) and coefficients L, J given by 


w[s^) — 


{£ + 1 ) ris ’ 


w{t‘) = 


l 8 (£ + 2 )(£ + 4)(£ + 0 )z(£) 


(£ + 3) 


’ 


The field equation for the reduced field is then 

" ^ ^ (A.37) 






where 


^^1^2^3 T^^1^2^3 ^^s^2 T ^^^3 T^i^ 2^3 g T^1^2^3 ^5^2 ^5^3 ‘ (A.38) 
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3^/^ 
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15^/^ 
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X4/22 


0 
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9\/T5 

0 

0 


0 

0 

0 

0 

5184 [2 

5 V 5 


Table 1: Coefficients {£>^ 22 ,-£^^>22?-^^> 22 }, {Ji> 22 ,L^ 22 } and A ^22 for the several 'tpg = {s^,t^} or 
'^i = {S^,T^} that appear in ( |A.39 )- |A..41| ). (Note that these coefficients depend on 
the harmonic representation we use for the S^). 


In practice, for our lumpy AdSsxS^ BH and up to the relevant order 0{z^), only the 
sources associated to contribute in the RHS of equations (A.34), ( |A.35| ) and (A.37). 

Accordingly, up to this order we can rewrite them simply as 

(□ - mlJ ^Pe = D^ 22 {s‘^? + E^ 22 D^s^D^^s^ + (A.39) 

4'^ = {ipe - J^ 22 {s^f - L^ 22 D^fD^f ), (A.40) 

= A^22(s^)^ (A.41) 


where the coefficients {D^ 22 ^ F'ij; 22 }^ {J'ip 22 ^ L^ 22 } smd A^22 foi* fhe several are given in 

Table 00 Explicitly, the 5-dimensional scalar fields of the lumpy AdS 5 x BH are 

5'=0, =-^y|/32!/+(^z^ + i(3 + (32!/pz‘‘) , 

=4 vf >'^3 =vf 4« ^ 

97 

7-0 = 0 , r' = 0, T^ = 0, T^ = 0, T^ = —^ylplz\ (A.42) 

A similar treatment for the massless KK graviton field, including a field redefinition, allows 
us to write the d = 5 graviton G/^iy in terms of the d = 10 fields. Again up to the relevant order 
{Xjz^')0[z^\ in the lumpy AdSsxS^ BH only contributions sourced by contribute. 

The reduced metric then reads (recall that is the AdSs nietric) 


^ __ 




9 


9 


9. 




This yields the following non-vanishing components 


1 


Gzz = ^ 


Gij — rjij 


65 

' 4^ 

O ~fo 




(A.43) 


(A.44) 


^ yX [85/3f K 2 + 30^2 Ki + 72Ki{lQ5o + S 4 - 192)] - 13824 ATi yl + 1152 
^ 18432 ^ 


“Note that these coefficients depend on the harmonic representation we use for the S° 


28 













































where rjij = diag{ — 1 , with being the line element of a unit radius and, to shorten 
the presentation, we introduced the auxiliary constants {Kq^ i^i, i^ 2 } such that {Kq^ i^i, i^ 2 } = 
{1,1,1} for i = j = t, while {Kq^ K 2 } = { — 1, —1/3, 31/17} for components i = j on the S^. 


The field equations (A.41) can be obtained from a 5-dimensional action, namely 


7Y2 r 

S5d=^j 






(A.45) 


with G/^iy given by (A.43) (recall that Newton’s constant G 5 is given by ( |3.2| )). The first two 
contributions in this action are the Einstein and cosmological terms (recall L = 1) that admit 
AdSs as a solution. Up to order O(z^), the potentials in this action are 




for = 5^ , 


- 48(r)^ = 0 {z^), for = 


otherwise. (A.46) 

Variation of the 5-dimensional action w.r.t. the scalar fields indeed yields the massive Klein- 


V^e) = 2"^* yy = 0{Z^), 


Gordon equations (A.41); note that □ = up to order 0{z'^). 

The Einstein equation that follows from the 5-dimensional action ( |A.45 ) is 


R/iA^] — 2 ( -h — -G/^iyT^ ) , 


(A.47) 


where is the Ricci tensor of and the energy-momentum tensor reads 




^ + vyi) 






(^^idsy + y(52)) + o(/). 


(A.48) 


As indicated by the second equality up to the relevant 0{z^) only the scalar field ^ con¬ 
tributes to the stress tensor. 


A.5 Holographic renormalisation and Stress tensor. 

We can now apply the standard holographic renormalisation procedure to the 5-dimensional 
Introduce the Fefferman-Graham coordinate Z — z — 


solution 


51 


65 

36864 


_ v\P 2 ^^^ 0 [z^) 

for the 5-dimensional metric and denote the boundary coordinates collectively by X. 

Furthermore, collectively denote the scalar fields of the system by $ = (where 

massive KK gravitons described at linear order by ( A. 20 )-( pV. 2 lD ) and recall that 
the conformal dimensions A (and A_ = 4 — A) of the operators dual to these fields that are 
given in Table 

The expansion around the holographic boundary Z = 0 for the 5-dimensional metric G^^y, 
and scalar fields $ = is then 


dZ‘^ 


— _ 

^2 ^ ^2 l^ij 


Gf;\x) + zyf^{x) + [Gf^{X) + log Z^H\f{X)) 


+ 


Z) = Z2 (log zHf^){X) + $2^) {X)+ ■■■), for A = Abf = 2; 
Z) = Z(4-^)$fo)(X) + •.. + Z^$f 2 ^_ 4 )(X) + •.. , for A > 2 , 


dX^dXR, 


(A.49) 
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5^ 

F 

(t>KK 

A+ = A 


£ + 8 

£ + 4 

A_ = 4-A 


-£-4 

-i 


Table 2: Conformal dimensions A± of the scalar fields and of the KK gravitons 


where the A = 2 case saturates the 5-dimensional Breitenlohner-Freedman (BF) bound, i.e. 
A+ = A_ = In Fig- 13, we plot the conformal dimensions of the several scalar fields 


as a function of the harmonic quantum number 1. The BF bound is saturated 




only for the field S 

In the above off-boundary expansion, the non-normalisable modes G(o)zj 5 ^(o)’ source 

terms for the boundary QFT stress tensor and dual operators of dimension A = 2 and A, respec¬ 
tively. On the other hand, the normalizable modes, namely determined 


■(0)’ ^(2A-4) 


by solving the field equations (A.47) and (A.41) of (A.45) subject to regular (ingoing) boundary 
conditions at the BH horizon. All other coefficients of (A.49), typically represented 


■{k) 


by dots, are expressed as a function of the (non-)normalisable modes and their derivatives. 





Figure 13: Conformal dimensions of the several operators that are present in the lumpy system. 

From top to bottom the points connected by a continuous line are the conformal 
dimensions A+ of the dual operator of the scalars T^, KK gravitons and 

scalars S^. The dots connected by a dashed line represent the conformal dimension 
A_ = 4 —A+; in the region displayed, the plot shows only the dashed line associated 
with the scalars S^. 

At this point we can discuss the boundary conditions (BCs) that we impose in the holo¬ 
graphic boundary Z = 0. We do not want to deform the boundary background so we fix 
to be the static Rt x metric as a Dirichlet boundary condition. Moreover, we do not want to 
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deform the boundary CFT, i.e. we require vanishing sources in the scalar fields, so we impose 
vanishing Dirichlet boundary conditions on the fields $. Altogether we thus have the asymptotic 
boundary conditions: 


g!-A) = 


^( 0 )W= 0 , 

^fo)W=0, 


( 0)1 

ij iRtxS^ 
for 

for 


A = Abf = 2; 
A > 2 . 


(A.50) 


These BCs can be discussed in more detail with the aid of Fig. IB A priori^ we could impose BCs 
that would allow operators with A > Aunit^ where Aunit = 1 is the unitarity bound. However, 
the dual CFT of our system is A/" = 4 SU(N) SYM. In this special case, the requirement that the 
norm of the supercurrent of the theory is positive definite requires that we exclude conformal 
dimensions in the range Aunit ^ A < Abf (ultimately this means that we have a dual SU{N) 
gauge theory and not U{N)] for further details see e.g. section 2.6 of [^). Therefore, we impose 
BCs such that only dual operators with conformal dimension equal or higher than the BF 
bound, A > Abf = are present. This means that the BCs we impose are such that sources 
of operators in (and below) the red area of Fig. [T^ vanish. 

In addition, some modes above this region (the source terms with {A_,£} = { 2 , 2 }, and 
{A_,^} = {4,0}) are eliminated by the BCs = 0 and = 0 on the scalars S^. There is 
however a special mode that we cannot exclude with our BCs namely, the mode with {A_, £} = 
{3,1}. This is a pure gauge mode and therefore it does not appear in any physical quantity]^ 
We cannot remove this mode since we are using the deTurck method. The gauge is fixed after 
solving the equations, and cannot be imposed from the equations of motion alone (see discussion 


associated to ( 2 . 1 )). 

We can now discuss the undetermined constants in the asymptotic expansion ( |A. 10 | ) that 


are permitted by the BCs (A.50). The coefficient in (A.49), with {A+,£} = { 2 , 2 }, is 


proportional to the parameter /32 appearing in ( |A. 10 | ). In ( |A.49 ), the normalisable modes 


^^A- 4 ) scalars associated to {A+,£} = {3,3} and {A+,£} = {4,4} are, respectively. 


proportional to the parameters 73 and ^4 present in ( A.10[ ). Finally, the fourth constant 5o 
that appears in the boundary expansion (A.10) describes a KK graviton with {A+,^} = {4,0}. 
Normalisable modes with A+ > 4 — see Fig. [T3]- appear only at an order in 2: higher than 


the one displayed in (A. 10). They are also present in our solution but we do not discuss them 
further because they do not contribute to the mass of the lumpy BHs. 

The normalisable modes are related to the holographic 1-point functions that give the VEVs 
(T^j), (O^) and {O^) of the dual operators, via the standard holographic renormalisation pro¬ 
cedure. In particular, the vev of the operators 0^2 and are (^(q) ^( 2 ) directly 


from ( A.42[ )) 


(G 53 ) 


/5 

^2*5(0) 7r2 8V3^+^^’ 

iV2 1 /5 o 

- ^^(2) - ^ 8V6^+^^’ 


(A.51) 


^Consequently, it is associated to a parameter that does not appear in the harmonic coefficients (A. 10) up to the 


order in 2 : that we display because (A. 10) shows only the order needed to compute the relevant physical charges 
and vevs. 
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and the holographic stress tensor is 


iV2 


+ 5'S(o)G!"> + ^ (ti g?,, - (1VOp|)=) G<; 


( 0 ) 

d 


(A.52) 


2^ 4 

which, for the lumpy BHs, explicitly reads 


^(G'?2))d + TrG(2) + 


iV2 


+ 4 y+ “ 2QJ2 ( ^^'^2 + 5/^2 + 12(16(5o + ^4 — 192) 


diag <; 1, - r?- 


(A.53) 


with r]--j being the line element of a unit radius S^. This holographic stress tensor is conserved, 
(o)v.(rd) = 0, and it is traceless, (T/) = oQ 

An important holographic quantity that we want to extract from (A.53) is the energy of the 


solution. This is done by pulling-back (Tij) to a 3-dimensional spatial hypersurface with 
unit normal n and induced metric -h , and contracting it with the Killing vector 

^ = dt that generates time translations. The integral of this quantity gives the desired energy 


E = - [ 
jjit 


Tin 


3072 


576 + 2304^(1 + ?/;)- 2/+ 5/^2 + 30/3^ + 12 (16 5o + ^ 4 ) 


(A.54) 


This is the main result of this Appendix. In the main text we will use (A.54) (rewritten in 


(3.4) with factors of L restored) to determine the energy of the lumpy BHs. The energy of the 
AdSs-Schw BH, EsAdS^/^^^ — (3/4)^^ (?/^ -h l) + 3/16, is recovered when we set the lumpy 
parameters to zero, {32 — 5^) — 5^ — 0. 


B Numerical details and validity 

In this appendix we discuss the validity of our numerical results, while giving further details of 
the numerical construction of the AdSsxS^ lumpy BHs. 

We start by testing the numerical convergence. We use pseudospectral collocation meth¬ 
ods, and thus we expect exponential convergence with increasing number of grid points. We 
demonstrate this convergence with the panel of Fig. In this figure, as a typical example of 
our results, we consider a lumpy BH at constant temperature T = 0.50065, and show how its 
entropy changes as the number N of grid points is varied. 

Next, we test numerical convergence of the norm of the deTurck vector defined below 
. The Einstein-de Turck method solves Einstein equations in the gauge — 0. Therefore, 
the norm of the deTurck vector is a measure of how well this gauge condition is satisfied, and 
verifies that we have a proper solution to the Einstein equations and not a DeTurck soliton 
with ^ 0. On the right panel of Fig. we take a lumpy BH at constant temperature and 
plot the square root of the norm of the DeTurck vector (evaluated at the asymptotic boundary, 
^ = 1, and at the rotation axis x = 1) as a function of the grid points. Again we confirm the 
presence of exponential convergence. We find that < 10“^^ everywhere. 

^The lumpy BH asymptotes to global AdSs which is conformal to the Einstein Static universe Rt x ^ and 
thus conformal to flat space. Therefore the gravitational conformal anomaly vanishes. Moreover a possible 
contribution, both to the conservation equation and trace, of the form is not present because we impose 

Dirichlet boundary conditions in the scalar field, = 0. 
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Figure 14: Left panel: Convergence test for the entropy of the lumpy BHs. We plot 
1 — |S'(A^)|/|S'(A^ + 2)1, as a function of the number of grid points N. Right panel: 
Convergence test of the square root of the norm of the DeTurck vector at the asymp¬ 
totic boundary and at the rotation axis, as a function of the grid points 

N. Both plots are for lumpy BHs with T = 0.50065 (i.e. = 0.44225). 


The first law and the energy read from KK holography provide a final important test of 
our numerics. We can independently compute the energy, entropy and temperature using KK 
holography. Therefore, we can test whether the first law, dE — TdS^ is satisfied. We find this 
to be the case as mentioned below (3.4) and in the discussions associated with Figs. |^and[^ 
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